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Abstract 

Geant4 low energy extensions have been used to simulate the X-ray spectra of industrial X-ray tubes 
with filters for removing the uncertain low energy part of the spectrum in a controlled way. The results 
are compared with precisely measured X-ray spectra using a silicon drift detector. Furthermore, this 
paper shows how the different dose rates in silicon and silicon dioxide layers of an electronic device can 
be deduced from the simulations. 



1 Introduction 

Radiation hardness of semiconductors can be tested 
by irradiation with X-rays, since these introduce ef- 
ficiently trapped ionization in the silicon oxide lay- 
ers, which e.g. leads to voltage shifts in MOSFET 
transistors. The damage depends not only on the 
X-ray flux, but also strongly on the spectrum of 
the X-rays. The total flux can be obtained from 
a calibrated pin diode, but an accurate measure- 
ment of the X-ray spectrum is complex and time 
consuming. A simulation is a more convenient way 
to get the spectrum. In this paper Geant4 [l] low 
energy extensions have been tested for this purpose. 
The results are compared with precisely measured 
X-ray spectra using a silicon drift detector (SDD). 
Furthermore, the different dose rates in silicon and 
silicon dioxide layers of an electronic device can be 
deduced from the simulations. 

2 Experimental Setup 

The X-ray tube of type DX-W8x0.4-2S with a tung- 
sten anode has a maximum power of 2000 W and 
a maximum voltage of 60 kV. The tube window is 
made of 0.4 mm Beryllium, so that only low en- 
ergy photons are absorbed by the window. The 
X-ray tube is equipped with different filters that 
can be mounted behind the exit window. These 
filters are vanadium (15 fjxa), manganese (25 /im), 
iron (15 /jm), nickel (15 /im) and zirconium (75 /im). 
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These filters absorb low energy photons in different 
energy regions. The minimal configurable current 
of the HV-module is 2 mA. The top part of Fig. [T] 
shows the tube with the filters. The tungsten colli- 
mator for measuring the spectrum at a reduced rate 
in a silicon drift detector is shown in the lower part. 

3 Dosimetry 

One critical information in the test for radiation 
hardness of semiconductor devices is the dose. This 
chapter describes the way of calculating the dose to 
any kind of material if the X-ray spectrum and a 
reference measurement are known. 

3.1 Dose rate calculated from X-ray 
spectrum. 

The X-ray energy spectrum provided by the source 
is given as N(E) which can be either measured 
or simulated. Since obtaining a flux normalized 
spectrum is even more difficult, N(E) is arbitrary 
scaled. Only a fraction of the spectrum is absorbed 
during irradiation. The absorbed part is deter- 
mined by: 

N abs (E) = N FUter (E) ■ (1 - e ^ E > d ), (1) 

where p represents the density, a the absorption 
coefficient and d the thickness of the material. The 
absorbed power can be calculated to: 

Pabs =a- / N abs (E) ■ E ■ dE, (2) 
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Figure 1: Sketch of the experimental setup. The 
lead collimators are used to mount the tungsten col- 
limators and reduce background radiation. 



where a represents the unknown scaling factor. The 
dose rate is defined as: 



D = 



Pa 



b.« 



(3) 



3.2 



Absorbed dose in a depleted 
diode as reference 



The missing parameter a can be obtained with a 
reference measurement. With a depleted silicon 
diode it is possible to measure the photo current 
from ionizing radiation. Each generated electron- 
hole pair represents 3.65 eV [2], so that the ab- 
sorbed power in the device can be calculated as: 



P a6ssi = ^*3.65ey, 



where e denotes the elementary charge. 
The dose rate in this device is then given by 



(4) 



D Sl = 



Pabssi 

msi 



(5) 



where msi is the mass of the diode. The measured 
dose rate is only valid for silicon material of the 
same thickness. In our measurements we used a 
planar diode with a thickness of 300 /im and a lead 
collimator with a circular opening of = 1.5 mm. 



Depletion over the whole thickness was ensured via 
a CV-measurement. 

Using Eqs.[2]and|4] one can now obtain a: 

Iphoto ■ 3.65eF 

a= , (6) 

e- J~N absSi (E)-E-dE 

where N abssi is the absorbed spectrum in the ref- 
erence diode. It is now possible to calculate the 
dose for any other material if the energy spectrum 
of the X-ray source is known. The spectrum can 
cither be measured or simulated. Both ways and a 
comparison are described in Sect. [4] 

3.3 Calculation of the dose rate in a 
different material 

Silicon dioxide is often used as an insulator in semi- 
conductor devices, which suffers from charge build- 
up caused by ionizing radiation. Hence, for radi- 
ation damage studies it is important to know not 
only the total dose rate, but also the dose rate in 
different layers, which usually have different absorp- 
tion coefficients. We use S1O2 as an example to 
show the calculation of the dose rate in any mate- 
rial. 

The absorbed power is calculated as in Eq. [2j and 
using Eq. [6] one obtains: 
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(7) 

The ratio of the absorbed dose in the reference ma- 
terial and in the target material is therefore used to 
scale the reference measurement P a bs Si to the ab- 
sorbed power in the target material P absS io ■ The 
dose rate in SiC>2 follows from Eqs. [3]and[iJ 

PabssiOo 



D s 



iO- 2 



m S i0 2 

Pabssi 

msi0 2 



J Ema * N absSiQ2 {E) ■ E ■ dE 
Jo ma * N abssi {E) ■ E ■ dE 



(8) 
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2y Cu 


47 Ag 


5e Ba 


«b Tb 


241 Am 




E[keV] 




E[kcV] 




E[kcV] 




E[keV] 


E[kcV] 


100 


8.048 


100 


22.162 


100 


32.193 


100 


44.481 


59.5 


52 


8.028 


53 


21.990 


54 


31.817 


56 


43.744 


13.9 


17 


8.905 


16 


24.942 


18 


36.378 


20 


50.382 








4 


25.456 


6 


37.257 


7 


51.698 








9 


24.911 


10 


36.304 


10 


50.229 





Table 1: List of materials used for the energy calibration an their characteristic X-ray energies. I re ; is 
the relative intensity of the peaks for the given material. 
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Figure 2: Energy deposits from the lines of the ma- 
terials in Table [TJ 
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Figure 3: Linear fit to the known energies of the 
observed lines in Fig. [2j 



The number of absorbed photons, N a j) SSi and 
N a bs S io 2 ' can ^ e calculated from Eq. [T] using the 
known energy dependent absorption coefficients a 
for each material [3]. In principle the method can 
be used for arbitrary materials and has been applied 
for all six possible filter settings. 

4 X-ray spectrum 

A measurement of the X-ray spectrum is often not 
possible, since the count rates are usually too high 
for detectors with a suitable energy resolution. Al- 
ternatively, the spectrum can be obtained from a 
simulation. In this chapter we describe both ways, 
simulations are done with Geant4, the measure- 
ments are done with a silicon drift detector. 

4.1 Measurement of the spectrum 

A silicon drift detector (SDD) with an integrated 
preamplifier of the type AXAS-P (KETEK GmbH, 



Munich, Germany) was used as detector. The sig- 
nals were further amplified by a shaping amplifier 
and registered by a readout system Multiport II 
(Canberra Industries, Inc., Meriden, CT, USA) 
with 16384 ADC channels. All measurements were 
done with an X-ray energy of 60 kV and 2 mA an- 
ode current. With this configuration the photon 
flux in the detector was too high to detect single 
photon events, since the flux was several order of 
magnitude above the count rate for the slow drift 
detector. Two tungsten collimators were used to re- 
duce the count rate. These collimators were 1 mm 
thick with a hole of ~0.1 mm diameter. 

One collimator was mounted at the exit window 
of the tube and another one on the detector, as 
shown in Fig. [T] 

A radioactive 241 Am source was used for the en- 
ergy calibration. This source includes different ad- 
ditional materials, which are excited by the 241 Am 
source and emit characteristic X-rays. The materi- 
als and their X-ray energies can be found in Table[l] 
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Figure 4: Geometry for the first simulation step, 
where all main parts of the X-ray tube are imple- 
mented. An electron beam hits a tungsten anode 
and the resulting photon spectrum is registered by 
a sensitive air layer (detector). 
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Figure 5: Geometry for the second simulation step. 
Photons were generated accordingly to the spec- 
trum generated by the tube configuration in Fig. [4] 
These photons are registered in a detector after 
passing a filter and a thin tungsten foil represent- 
ing the absorption by scattering on the inside of the 
collimator in Fig. [T] 



The relation for the energy calibration is: 

E[eV] = (5.04 ± 0.003) • ADC + (40 ± 19), 

as obtained from a linear fit to the observed energy 
peaks in ADC channels, shown in Fig. [3] The small 
errors show that the calibration is precise. 

The SDD has a thickness of 300/im. Some pho- 
tons pass through the detector, especially at higher 
energies. To compensate for this effect the absorbed 
spectrum was corrected by weighting it with the in- 
verse photon absorption probability (3| for a layer 
of 300 fim silicon. The energy resolution of the 
measurement was determined by fitting a Gaussian 
curve to the L a peak. The full width at half max- 
imum (FWHM) is 0.5 keV, which is comparable to 
the resolution 0.58 keV obtained during the energy 
calibration. 



Figure 6: The simulated spectrum of the tungsten 
anode using the setup schematically presented in 
Fig. [4] (upper curve) and the reduction by the thin 
tungsten foil representing the absorption by the col- 
limator in Fig. [5] without filter (lower curve) . This 
lower curve is compared with data in Fig. |7(a) The 
La, L/3 and L 7 lines (8.34 keV, 8.4 keV, 9.67 keV, 
9.96 keV, 11.3 keV) of tungsten as well as the con- 
tinuum of bremsstrahlung are clearly visible. 



4.2 Geant4 simulation 

Geant4 is a toolkit developed at CERN for the simu- 
lation of the passage of particles through matter [l] . 
It is an object-oriented simulation framework which 
provides a diverse set of software components [4]. 
All aspects of the simulation process like the ge- 
ometry of the system, tracking of particles through 
materials and external electromagnetic fields as well 
as the response of sensitive detector components are 
included in the simulation. Geant4 provides also a 
set of different physics models to describe the inter- 
actions of particles with matter across a wide energy 
range. The standard electromagnetic physics pack- 
age provides implementations of electron, positron, 
photon and charged hadron interactions. It is suit- 
able for most of the Geant4 applications, but it is 
not optimized for low energy particles. Two specific 
low energy electromagnetic models are available: 
Livermore (based on the Livermore library [5jj7]) 
and Penelope [8J. Geant4 version 9.3.p01 was used 
and the version of CLHEP libraries was 2.0.4.5. At 
first both low energy physics packages were used. 
However, results obtained with Livermore showed 
differences in comparison with measured spectra. 
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The characteristic tungsten lines had a lower inten- 
sity compared to the bremsstrahlung slope. The 
results presented here arc obtained with the Pene- 
lope package. 

In order to reduce the simulation time the pro- 
gram was split into two different parts. In the first 
part the pure tungsten spectrum of the X-ray tube 
was generated. Then this spectrum was used as 
input for the second part, where the effect of dif- 
ferent filters and collimators could be studied. The 
first step of the simulation is the generation of an 
electron beam with an energy of 60 ± IkeV. The ge- 
ometry is shown in Fig. [4j The environment is vac- 
uum. The vertical position of the beam was varied 
by 12 mm to simulate the glow filament (the elec- 
tron source within the X-ray tube). The electrons 
were shot at a tungsten anode at a distance of 8.5 
mm. Below the anode a 0.4 mm thick beryllium 
window was placed. After the photons, generated 
at the anode, passed through the beryllium window, 
their energy was registered by a sensitive air layer 
(detector). 

In the following step this photon spectrum was 
used as input for the setup shown in Fig. [5] Mate- 
rial and thickness of the filter were taken from the 
X-ray tube specifications. After the filter a tung- 
sten collimator was used to collimate the beam. 
The hole geometry of this collimator is irregular, 
since the small hole of ~0.1 mm was produced by a 
laser shot. The scattering on the inner edges of the 
two tungsten collimators was simulated by a thin 
tungsten foil of 3 /zm placed between the filter and 
the detector. The thickness was chosen to agree 
with data, as will be shown in the next section. In 
Fig. [6] the simulated spectra before and after this 
thin tungsten layer are shown. 

4.3 Comparison of Measurements 
and Simulations 



Filter 


Dose rate 


Dose rate 




simulation 


RadFETs 




Gy/s 


Gy/s 


Zr 


0.24 


0.22 


V 


1.84 


1.93 


Mn 


0.86 


0.90 


Fc 


1.32 


1.19 


Ni 


1.44 


0.99 


Air 


6.29 


5.95 



Table 2: A comparison of the simulated and mea- 
sured dose rates after several filters in RadFets. The 
values correspond to a voltage of 60 kV, 33 mA tube 
current, a distance between tube and the object of 
155 mm and a SiC>2 thickness of 130 nm. 



expected, since in the simulation photons are regis- 
tered, not detected by the detector. Therefore, the 
simulation does not include additional absorption 
effects in the detector, e.g. by Compton scattering. 
The small enhancement at 25.2 keV, which appears 
in every measured spectrum, is the K a line of tin [9] . 
Characteristic X-ray photons from the solder used 
to connect the peltier cooling in the detector cap- 
sule are responsible for this peak. 



5 Dose rate measurements 
with RadFETs 

To make sure that the results are reliable we cross- 
checked the simulation with special ionizing radia- 
tion detectors, namely RadFETs, which consist of 
calibrated MOSFETs, where the shift in the thresh- 
old voltage by the positive charge build-up in the 
silicon oxide layer has been calibrated 



10 11 . Ta- 



For a realistic comparison the simulations were 
convoluted with a Gaussian curve with the same 
FWHM as the detector. The measured energy spec- 
tra after all available filters are in excellent agree- 
ment with the simulations, as shown in Fig. [7] 
Note the reasonable agreement between simulation 

and measurements without filter (left top panel), 6 Summary 
which is the justification for the thin tungsten foil in 
Fig. [^representing the absorption in the collimator. 
Only the low energy part differs slightly, which is 



ble [2] compares the results from the simulation with 
the measurements. One observes that filtering out 
the soft part of the spectrum (corresponding to the 
lower dose rates) improves the agreement between 
measurements and simulations, as expected since 
the low energy part of the spectrum in less well 
simulated, as discussed before. 



This paper shows a convenient method to obtain the 
applied dose from X-ray irradiation in the various 
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Figure 7: A comparison of the simulated and measured spectra with different filters. 
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layers of an electronic device, namely by one calibra- 
tion measurement of the total flux and a simulation 
of the X-ray spectra, from which the dose rates in 
materials with different absorption coefficients can 
be calculated. The method has been crosschecked 
with calibrated RadFets. The difficult part is the 
reliable simulation of the X-ray spectrum after the 
usual filtering. We found that the Geant4 X-ray 
simulation with the Penelope physics package is in 
excellent agreement with precise measurements of 
the spectra using a silicon drift detector. 
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